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ABSTRACT 


-We— consider  a  system  consisting  of  two  modules,  one  of  which  cannot  be 
easily  inspected,  while  the  other  is  monitored  continuously.  A  fault  devel¬ 
oping  in  the  first  model  is  called /’unrevealed"  (U)  while  the  one  in  the 
second  module  (which  is  assumed  to  be  detected  immediately)  is  calletj/^fe- 
vealed*'  (R).  It  is  supposed  that  repairs  are  initiated  as  soon  as  an  R 
is  observed,  but  not  otherwise.  Phi-1 14  p»  (Mi crdol e ct-ron-r-Re liatr.  (19 79 ) 
4J&r^4P5-504^--Reiriab7--gagineer-  (4S1&1) ,  ^221-231)  assumed  that  both  modules 

are  always  replaced.  -In^Ihis  paper  w*r  allows for  replacement  of  the  first 
module  only  if  a  U  fault  is  found  on  special  inspection.  Expressions  are 
derived  for  the  joint  distribution  of  time  to  first  repair  and  time  between 

■'i 

first  and  second  repairs.  Special  cases  aTe  considered  in  which  fairly 
simple  results  are  obtained  by  general  arguments  .without  recourse  to  anal¬ 
ysis.  Detailed  formulas  are  developed  for  a  particular  parametric  model. 
Relationships  with  relevation  theory  are  indicated. 
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□  □ 


1.  Introduction 


Phillips  [1,2]  has  considered  a  system  consisting  of  two  modules,  one 
of  which  cannot  be  easily  inspected,  while  the  other  is  monitored  contin¬ 
uously.  A  fault  developing  in  the  first  module  is  called  "unrevealved  (U)  - 
that  is,  until  a  special  inspection  is  carried  out  -  while  a  fault  in  the 
second  module,  which  it  is  assumed  will  be  detected  immediately,  is  called 
"revealed"  (R) .  It  is  supposed  that  repairs  are  initiated  as  soon  as  an 
R  is  observed  but  not  otherwise.  As  in  [2]  we  will  further  assume  that  re¬ 
pairs  are  effected  instantaneously  and,  if  carried  out,  result  in  the  re¬ 
paired  module  being  "as  good  as  new". 

We  will  be  especially  concerned  with  the  distribution  of  the  time  at 
which  a  second  (or  later)  repair  is  needed,  on  similar  assumptions  to  those 
of  Phillips,  except  that  we  will  allow  for  the  possibility  that  repairs  are 
made  only  to  modules  with  faults.  This  would  mean  that  if  an  R  occurs,  the 
first  module  is  specially  inspected,  and  replaced  only  if  it  has  a  U  fault. 

Phillips'  model  is  based  on  three  random  variables:  X:  time  from  repair 
of  R  to  next  occurrence  of  R,  assuming  no  U  present.  Y:  time  from  repair 
of  U  to  next  occurrence  of  U,  assuming  no  R  present.  Z:  time  from  a  U  fault 
to  an  R  fault.  X  and  Y  are  assumed  independent.  Z  is  independent  of  X  in 
both  [1]  and  [2],  but  may  depend  on  Y  in  [1],  though  not  in  [2].  We  will 
retain  this  possible  dependence. 

Using  the  notation  f^(w)  to  denote  the  density  function  of  a  random 
variable  W,  and 

00 

V*)  «  /  fw(t)dt 

to  denote  its  survival  function,  the  density  function  of  the  first  replace- 
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ment  time,  Tj,  is  ([1],  equation  (1)) 

1 1 

f“T  (tj)  -  fx(tl^SY(tl5  +  /  Sx(t1-2)fy>2(t1-z,z)dz  (1) 

(We  use  subscripts  to  denote  the  order  of  the  replacement  period.) 


2.  Distribution  of  time  between  first  and  second  replacements 

If  "repair"  consists  of  complete  replacement,  as  in  [1]  and  [2],  then 
the  time  from  first  to  second  replacement  has  the  same  distribution  as 
Tj,  and  and  are  mutually  independent.  Hence  the  time  of  second  re¬ 
placement  (Tj+Tj)  is  distributed  simply  as  the  convolution  of  two  Tj's.  Ex¬ 
tension  to  later  replacement  times  is  straightforward. 

The  situation  is  different,  however,  if  the  first  module  is  replaced 
only  if  a  U  fault  is  found  to  be  present  on  inspection  after  an  R  fault 
occurs.  When  such  a  fault  is  not  present  we  start  the  second  replacement 
period  with  a  first  module  already  aged  Tj,  rather  than  a  new  one.  Ke  note 
in  passing  that  the  same  situation  arises  when  replacement  of  both  modules 
is  automatic,  when  an  R  fault  occurs,  if  the  first  module  is  replaced  from 
aging  stock,  with  aging  occurring  at  the  same  rate  in  storage  as  in  service. 
Relevation  theory,  which  we  have  discussed  in  [3]  and  elsewhere,  is  appli¬ 
cable  to  this  situation. 

If  a  U  fault  is  present,  then  X2>  Y2  and  Z2  will  have  the  same  joint 
distribution  as  XJf  Y  and  Z^  and  the  sets  of  variables  (Xj.Yj.2j)  and 
(X2,Y2,Z2)  will  be  mutually  independent.  If  a  U  fault  is  not  present,  then 
given  Tj  ■  tj,  we  will  have 

X2  distributed  as  -  density  function  fx(x2l; 

Y2  distributed  as  conditional  on  Yj>tj:  -  density  function 

fY(y2+tl)/Sy(tl}  and  survival  function 
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Z2»  given  Y2>  distributed  as  Zj  given  Yj»Y2>  so  t*lat  the  )oint  density 
function  of  and  Z2  is 


fZ|yCz2|y2Kv(y2IV,i) 


(y2,22) 


(y2) 


VVlii 

MV 


(2) 


Also  we  have 


P^tp  *  Pr[U  fault  present (T^tj]  =  l-f^tjJSy (tj)/fT  (3) 

The  conditional  density  function  of  T2>  given  T^*^,  is  therefore 

■  PU(tl){fX(t2,SY(t2)*42sx(t2'2)fy,Zl,2'J’*,dzl 


'-pu(ti> 


fY(tl4t2"2-1 


+  Sy(t1)  {fX(  WW^O  SX(t2“Z)fY,Zft2"Z,Z)  fy(t2*z)  d=*  (0<t2) 

(4) 

The  joint  density  function  of  T^  and  T2  is 


“T1,T2(tl,t2)  =  fT2lT1<‘t2 fTj(V 

t 


{fX(t2)SY(t2)4/o2sX(t2-z)fY,Z(t2'Z'z)dz} 

*  ^0  SX^tl"Z^fY,Z^tl~Z,Z^dZ 

4  {fX(t2)SY(tl4t2)+^0  SX(t2"Z)fY,Z(t2"Z,Z)  fy(t2>z)  dZ^W 


fy(tj4t2_z) 


{fX(t2)SY(t2)4H(t2'°)}H(tr0) 


♦  {^WW  4  H(t2'tl)}fX(tl)  (°<tl‘t2) 
where  H(a,b)  *  Jj  Sx(a-z)fy  Z(a-z,z){fy(a+b-z)/fy (a-z) }dz 

The  overall  density  of  T2  is 

t, 

fT  (t2)  *  {fx(t2)Sy(t2)+/Q  Sx(t2-z)fy^z(t2-z,z)dz} 


(5) 

(5a) 

(ft) 


r0  /0lsx(tr2)fY.2(,rz>2)dId,i 
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+  fX(t2),^0fX(tl)SY(tl+t2)dtl+^0fX(t1^02sx(t2‘z)fY,2(t2-Z,Z)  Yfy(t2-z)  dz 


We  note  that  /^S^t^-zjfy  ^(tj-ZjZjdzdtj  is  simply  the  overall  proba¬ 
bility  that  a  U  fault  is  present  at  the  first  replacement  time,  that  is, 
Pr[X>Y] ,  while  f°Q  VVWV  =  Pr[Y>X-t2]. 


3.  Some  special  cases 


Before  proceeding  to  study  cases  in  which  the  joint  distribution  of 
X,  Y  and  Z  is  completely  specified,  we  first  take  note  of  a  few  results  of 
broader  character  which  can  be  derived  by  general  reasoning,  without  recourse 
to  analysis. 

Intuitively  one  might  feel  that  1  -  P^(tj)  =  SyCtj)*  since  the  proba¬ 
bility  of  no  U  fault  occurring  up  to  time  t,  is  Sy(tj),  but  this  does  not 
take  into  account  the  differential  effect  of  the  times  of  occurrence  of  R 
faults  according  as  they  are,  or  are  not  preceded  by  U  faults.  If  the 
time  to  an  R  fault  (T^X  if  XsY;  *  Y-»Z  if  Y<X)  were  independent  of  Y,  then 
f?  (tj)  s  fj,(tj)  and  so,  we  would  have  1  -  P^ftj)  *  Syftj).  It  is,  however, 
clear  that  this  is  not  usually  the  case;  indeed  independence  of  Z  and  Y 
does  not  imply  independence  of  Tj  and  Y;  in  fact  it  usually  implies  depen¬ 
dence. 

If  replacements  of  first  modules  are  from  aging  stocks,  with  aging  being 
the  same  whether  in  storage  or  in  service,  ("relevated"  using  the  terminology 
of  [37]),  then  it  makes  no  difference  (so  far  as  distributions  of  T\’s  are 
concerned)  whether  or  not  a  first  module  is  replaced  when  there  is  no  U  fault. 
This  is  because  after  the  repair  we  have,  in  either  case,  a  first  module  of 
age  tj;  and  this  result  is  valid  whether  replacement  of  second  modules  is 
from  aging  stocks  or  not. 
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If  the  lifetime  distribution  of  the  first  module  (that  is,  of  Y) 
is  exponential,  then  (X^Y^Zj)  and  (X2,Y2,22^  and  so  Tj  and  T 2  will  be 
mutually  independent,  whether  replacement  of  the  first  module  is  from  aging 
stocks  or  not  (or  is,  indeed,  performed  or  not  when  no  U  is  present).  This 
follows  from  the  lack  of  memory  (or  old-as-good-as-new)  property  (see 
e.g.  t4l)  the  exponential  distribution.  If  replacement  of  the  second 
module  is  from  aging  stock,  this  may  no  longer  be  the  case,  though  it  is 
so  if  the  distributions  of  X  and  Y  are  both  exponential,  whatever  the  dis¬ 
tribution  of  Z,  provided  this  depends  only  on  time  (Y)  since  replacement, 
and  not  on  actual  age. 

If  the  joint  distribution  of  Y  and  Z  is  of  the  Farlie-Gumbel-Morgenstem 
form,  with 

fY,z(X.z)  *  fy(y)f2(2)[l4a{2Sv(y)-l}{2S2(z)-l}]  (|o|<l)  (8) 

(see,  e.g.,  [5]),  then  in  (5a) 
a 

H(a,b)  =  /  Sx(a-z)fY(a*b-z)fz(z)[l+a{2Sy(a-z)-l}{2Sz(z)-l}]d= 

(9) 

When  we  come  to  consider  explicit  forms  for  the  distribution  functions, 
numerical  evaluation  of  the  distributions  of  Tj  and  T2  is  straightforward, 
but  for  elegant  analytical  results  there  are  technical  difficulties,  cen¬ 
tering  mainly  on  the  form  of  Sy(y)  and  in  particular  of  the  ratio 
fy(Vt2-z)/fy(t2-z),  and  also  on  the  form  of  Sz(z).  Assuming  exponen¬ 
tial  forms  does,  as  we  have  seen,  give  simple  results,  but  these  are  per¬ 
haps  too  simple  to  make  good  examples.  We  will  take,  as  an  illustrative 
special  case,  X,  Y  and  Z  to  have  density  functions 

ej2t  exp(-te*1)  (o<e.,t) 
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i«l,2,3  respectively,  with  the  joint  density  function  of  Y  and  Z  of  the 
Farlie-Gumbel-Morgenstero  form  (8) .  We  have 

fY(Wz)/fY(t2'z)  *  (tj+tj-zHtj-z)  e  ^;Sx(x)  «  (l+xe^Je  1 

-1  'y02  -1  ‘2/63 

Sy(y)  =  (l+y02  )e  *1  S2(z)  =  (l+203i)  e  3. 


Now  from  (5) 


-2  -1 
\  ‘■a  .  a  a 


fT  ,T  ^I’V  *  iei  t2^62  t2)  exp{-t2(61%62i)+H(t2,0)]  Hft^O) 

•t  /  0 

♦  [6j2f2{i+e‘1(t1+t2)}exp{-t2(e‘1+e21)-t1e‘1}+H(t2,t1)]e’2t1e  1  1  (io) 


with  H(a,b)  =  (62®3^  2  ®xp{-a(0j*«'02*)-be2*}/{l+(a-z)6“1}z(a+b-z)e 

a  .  -(a-2)/0  -2/0 

*  [l+a-2a{(lA-±)e  z+(z/0„)e  3) 

°2  •* 

-(a-2)e:1-z031 

♦  4a(l+=gj)  (z/dje  *  0  dz 


z(0i1+021-031 


(ID 


Although  (11)  is  complicated,  it  involves  only  integrals  of  form 
/a  zaexp(Bz)dz  *  B-a-1/^wae‘wdw  and  so  can  be  expressed  in  terms  of 
elementary  functions. 

Similar  analyses  can  be  performed  for  other  choices  of  distribution 
for  X,  Y  and  Z  (e.g.  Weibull)  but  they  lead  to  even  more  complicated  ex¬ 
pressions  than  (11).  We  hope  to  provide  some  detailed  analyses  of  such 
cases,  in  later  work,  and  also  to  introduce  cost  functions,  on  the  basis 
of  which  replacement  strategies  can  be  compared. 
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